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Abstract 

The scalar partner of the top quark (the stop) is relatively light in many models 
of supersymmetry breaking. We study the production of stops at the Large Hadron 
Collider (LHC) and their subsequent decays through baryon-number violating cou- 
plings such that the final state contains no leptons. A detailed analysis performed 
using detector level observables demonstrate that stop masses upto ~ 600 GeV may 
be explored at the LHC depending on the branching ratios for such decays and the 
integrated luminosity available. Extended to other analogous scenarios, the analysis 
will, generically, probe even larger masses. 

The status of global symmetries in particle physics is a much debated one. Of par- 
ticular relevance are fortuitous symmetries such as baryon number {B) and lepton number 
(L) . Apparently a consequence of the particular choice of the field assignments within the 
Standard Model (SM), each is left intact to all orders within perturbation theory, only for 
a combination to be broken by nonperturbative effects. While guaranteeing the stability 
of the proton, this leaves no room for the observed baryon asymmetry in the universe. As 
is well known, the latter needs, apart from enhanced levels of CP violation and a phase 
transition (or, at least, a non-equilibriated state of the universe), additional sources of 
S-violation (or L-violation for leptogenesis) . 

This has led to a sustained study of models of physics going beyond the SM ad- 
mitting 5-violation. Particularly well-motivated scenarios pertain to grand unification 
and supersymmetry (SUSY), whether considered separately or in conjunction. Indeed, 
low energy supersymmetry is considered to be, perhaps, the most attractive extension of 
the SM. The most general superpotential of the minimal supersymmetric standard model 
(MSSM) contains, apart from the generalizations of the Yukawa terms and the Higgs 
potential of the SM, additional terms violating B and L. Since, this would, nominally, 
lead to rapid proton decay, such terms need to be eliminated, and a global symmetry was 
introduced [T],|2] for this very purpose. Subsequently, a discrete symmetry, i?-parity, was 
shown to prohibit all dimension-four B- or L-violating terms |3]. However, as has long 
been realised, proton decay can be eliminated even when i?-parity is broken provided at 
least one of B and L are conserved. This is a particularly fascinating situation, for it 
admits a whole host of phenomenological consequences absent in the i?-parity conserving 
MSSM. Indeed, the collider signatures of supersymmetric particle production now under- 
goes a sea-change, for the lightest of the S'M-superpartners is no longer stable. And while 
it has been argued that the introduction of i?-parity violation would deprive the MSSM 
of one of its most atractive features, namely a natural Dark Matter candidate, note that 
the gravitino could yet be stable on cosmological scales. 

In this paper, we shall concentrate on the case where B is broken, but L is not. 
Such a situation is well motivated in a class of supersymmetric grand unified theories 



(GUTs) [4]. This, then, allows for terms in the superpotential of the form 

Wp^ 3 X';^,W^&^D'k , (1) 

where U}^ and Z)]j, denote the right-handed up-quark and down-quark superfields respec- 
tively and the couplings Aj^^ are antisymmetric under the exchange of the last two indices. 
The corresponding Lagrangian can then be written in terms of the component fields as 

= [n^dl + u^d]dl + u*d^^dl) + h.c, (2) 

thus allowing a squark to decay into a pair of anti-quark^ violating B and R. 

At a hadron collider, the couplings of eqn.(l2]) would lead to resonant production of 
a scalar, with the rates being potentially large if at least one of the quark fields belongs 
to the first generation. This could, in principle, lead to spectacular signatures both for 
dijets (invariant mass spectrum) or more complicated final states OE]- However, many 
of the A" couplings are severely constrained by a variety of low-energy observables [3-[9] , 
especially for operators involving at least two fields of the first generation|§. Given these 
two contrasting pulls, we restrict ourselves to the (A"-independent) pure supersymmetric- 
QCD processes as far as squark production is concerned, and consider a role for A" only 
in their decays. 

The mechanism of SUSY breaking is open to speculation. However, most such mech- 
anisms postulate that the sfermion masses be unified at some high scale, the extent of 
unification being mo del- dependent. This has the additional benefit of suppressing poten- 
tially large flavour changing neutral currents. In coming down to the electroweak scale, 
though, the masses would suffer substantial renormalization group evoulution. Further- 
more, the large Yukawa coupling of the top quark plays a role not only in this evolution, 
but also in engendering a mixing between the two super-partners of the top quark (the 
stops). Together, these effects very often result in the lighter stop being rendered the 
lightest of the strongly interacting superpartners. 

Thus, it is phenomenologically interesting to search for top squarks in the hadron 
colliders in the following scenario: (almost model independent) pair production of titi 
via strong interaction and decays thereafter via i?-parity violating process controlled by 
the parameter A |10| . 

It may be noted that there are scenarios beyond the SM other than SUSY in which 
a strongly interacting elementary particle may decay into a pair of quarks. GUTs, for 
example, abound in these [TT]. Such "diquarks", of which squarks in the i?-parity violating 
theory constitute but one example, can occur in various forms. They could be Lorentz 
scalars or vectors, SU{2)l singlets, doublets or triplets and, under, SU{3)c transform 
either as a 3 or a 6. While only some of the above combinations could appear in a trilinear 
coupling with two quark fields, others can participate if a quadrilinear term involving, say 
the Higgs, is considered jT2]. In the latter case, although the effective diquark-quark- 
quark vertex would exist only after the electrweak symmetry breaks, there is no a priori 
reason to relate the coupling to the quark mass terms. Diquarks, in particular, have 

^ Unless otherwise mentioned, for every process, the charge conjugated process is also implied. 
^Some of the couplings involving the third generation can be quite large though. 
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Figure 1: Production cross-sections for scalar (left) and vector (right) particles at the 
LHC for different collision energies. 



been shown [13] to be leading contenders for explaining the anomalously large forward- 
backward asymmetry in top production as reported by CDF |14| . 

Yet another class of strongly interacting particles may decay into a pair of jetfl 
namely color-octets like axigluons jT5], colorons [16] or Kaluza-Klein excitations of glu- 
ons [1^ or even color-triplets like excited quarks [T^. Such objects, though, can be 
formed as resonances, often with unsuppressed couplings, rendering the search for them 
to be much simpler [T91I20] . We will not consider such possibilities and would limit our- 
selves to the more difficult task of diquarks (squarks) with a small coupling to a pair of 
quarks. 

As is well-known, considerations of the p-parameter stipulates that the mass-splitting 
between members of a SU{2)l must be small. Thus, all other quantum numbers being 
equal, the SU(2)l singlet diquark would have the smallest QCD production cross sec- 
tion. Similarly, for all coloured SU{2)l singlet scalars, the one in the SU{3)c triplet 
representation would have the lowest production cross-section. 

As for vector diquarks, various inequivalent choices for the kinetic term (and, hence, 
for the coupling with the gluon) are possible (see Ref. [21]). Restricting ourselves to the 
case of minimal coupling, we find that the resultant cross-sections are significantly larger 
than is the case for the corresponding scalar (see FigH]). Thus, the choice of the squark 
as our template diquark is the most conservative one. 

Of course, the signal size at a collider detector depends not only on the total cross sec- 
tion, but on the kinematic distributions as well. It is worthwhile to note — see Fig|2U^a, b) — 
that the scalar and vector diquark have very similar angular distributions. In other words, 
the rapidity-dependence of the efficiencies would be very similar for the two cases. On 

^Distinguishing decays into a quark-antiquark pair or into a quark-gluon pair from a decay into a 
quark-pair (or an antiquark pair) is possible only if both the daughters carry a heavy flavour. 

''Note that while this choice results in a smaller production cross section than is the case for a Yang- 
Mills type coupling, it is possible to get marginally smaller cross sections by tuning the anomalous dipole 
and quadrupole moments. 
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Figure 2: Normalized kinematic distributions for QCD-driven pair-production of a diquark 
pair, (a) the rapidity; (6) the rapidity gap; (c) the transverse momentum; and [d) the 
invariant mass. 



the other hand, the spectrum is decidedly harder — see Fig|2]^c, d) — for the vector case, a 
reflection of the momentum dependence of its coupling to the gluon. Since we would need 
to impose strong p-r-cuts, the overall efficiency would be better for the vector case. Thus, 
on this count too, the choice of the squark as a template is a conservative one. 

Several searches — both in R-parity conserving and violating scenarios — have been 
made at LEP [22] and Tevatron |23] without success. These searches claim to exclude, de- 
pending on the relevant SUSY parameters, rrif upto ~ 200 GeV assuming 100% branching 
ratio for the respective decays. Such bounds, of course, weaken in the presence of other 
decay channels |24] . 

The high statistics data expected from the Large Hadron Collider (LHC) at -y/s = 
14 TeV provides an attractive option of looking for the ti in its i?-parity violating decay 
mode(s). At the partonic level, our signal events would read 

pp titi* {bq){hq) ■ 

We use PythiaG |2^ to generate both signal and background ( tt, W~^W~, ZZ ) events. 
Pure QCD background events (66+2-jets) have been generated at the parton level with 
Alpgen2 and interfaced with Pythia. 

In Pythia, partons in an event go through the quarks from hard processes undergo 
parton shower development and hadronization followed by decays. Pythia also emulates 
the initial (ISR) and final state radiation (FSR) from the partons for all the event samples. 
Tools provided by Pythia have been used to define a toy calorimeter with the broad fea- 
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tures of the LHC detectors: pseudo-rapidity coverage (—4.5 < rj < 4.5) and segmentation 
(A?7 X A0 = 0.1 X 0.087) [27]. To simulate the finite energy resolution, the jet momenta 
are smeared with a Gaussian function with a energy- dependent width 6E given by 



E ( GeV) 

with the two contributions being added in quadrature. The final state particles are passed 
through this toy calorimeter for forming jets. We also look for leptons {i = e, fi) but within 
the more restricted tracking coverage of the LHC detectors. Since leptons are measured 
quite accurately, we use the generated values of their energy, momentum and direction. 
Events are reconstructed in the following manner: 



• Calorimeter cells with Et > 1 GeV are used for jet reconstruction. The cone algo- 
rithm with, AR = A/(Ar/)2 + (A0)2 = 0.7, has been used requiring E^"* > 30 GeV. 

• Leptons with —2.4 < rj < 2.4 and Et > 5 GeV are deemed identified. 

• A jet originating from a 6-quark (&-jet) is identified by the decay length of the B- 
hadron(s) within it. The efficiency for tagging 6-jets is tuned to the expected value 
(~ 50% for 6-jets in tt events). 




For the signal events, we expect four jets with high E^ and no lepton with a high pt- 
We should not expect any except for that arising due to mis- measurement of the jets 
and neutrinos from secondary and tertiary decays of long lived particles. There may be 
additional objects (extra jets, I^t, leptons, etc) arising from underlying events (a serious 
concern at the high luminosity runs of the LHC) as well as due to ISR and FSR. But such 
objects are expected to fail hard cuts used in this analysis. At this juncture, we could have 
considered imposing an upper cut on I^x to remove, say, the tt background. However, 
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Figure 4: Distribution of Mjj and M^i, for signal events (rrif^ = 550 GeV) and two major 
background sources after the cut on the number of ordinary and b-tagged jets has been 
applied. 



since the jets from the ti decays have a high Ex, even a small relative mis-measurement 
could result in a substantial I^t leading to the rejection of a large fraction of the signal 
events (see Figj3]). Given this, and the relative unimportance of the ti background, we 
choose not to impose such a criterion. 

Instead, we use the following criteria for selecting signal events and rejecting back- 
grounds. 

1. There should be only four jets in the event {Njf^t = 4) with E^^ > E^^:^^. Of these, 
two are tagged b-jets 62) while the other two (ji, 72) are ordinary jets, i.e., they 
do not satisfy b-tagging criteria. (Sel 1). 

The jets are ordered descending in Et- 

2. Events containing any isolated lepton with pt > 15 GeV are rejected. (Sel 2) 

3. We also require ^bib2i'^jij2 ^ {"^^i 100} GeV ( see plots in FigH] ) to reject events 
with genuine W 01 Z (Sel 3). The plot for Mjj for tt events in Fig. Hlexpectedly 
shows the W mass peak. 

4. To eliminate the huge soft QCD background which may mimic signal with fake 6-jet 
tagging, we require (see right plots in Fig. [3]) ^ -E^^* > S'min- (Sel 4) 

For each value of m^- , nine different combinations of Sel 1 and Sel 4, viz, E^^^^ = 
{30, 50, 75} GeV and ^min = {400, 500, 600} GeV along with Sel 2 and Sel 3, were tried 
and the combination that gave the best signal visibility has been used for the final result. 
It may be noted that, in this analysis, we have used calorimetric observables which may 
be estimated fairly well both in magnitude and direction using parametric simulation. 

The events are sought to be reconstructed by computing the invariant masses of a 
pair of jets, one being tagged a -B-jet, and the other not so. Clearly, four such pairings 
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Figure 5: The left plot shows distributions ( normalized to arbitrary luminosity ) ofm'^gi^^ 
(black dashed line) and m'^^^^g (red solid line ) for signal events fm^~ = 550 GeVj along 
with m'^^gf^^ after Sel 5 (black solid line). The right plot shows distributions of m'^^gf^^ 
(normalized to shape) after Sel 5 for two major background sources. 



occur, viz. {biji,b2j2,bij2,b2ji}. Either 62^2 (combination 1) or, bij2,b2ji (combi- 
nation 2) represents the right combination in terms of the quarks from the hard processes 
(the decays) while the other represents the wrong combination. The two invariant masses 
belonging to the right combination would be expected to have a difference smaller than 
that for the wrong combination. To further reduce the chance of choosing the wrong 
combination due to accidental large mismeasurement of jet energies, we require (Sel 5) 

Am = mm{|mfejj^ — mb^j^l, Imb-^j^ — mbjjj} < 20 GeV . 

The average mass of the right combination (mb^j^ + mi,^j^)/2 or {mi,^j^ + mfe2jJ/2 

as specified by Sel 5 is then plotted in Fig|5l The other average is named m'^'J^^^g. The 
plot for m'^igfi^ in signal events shows the peak close to the input value of while the 
plot for m'^^g^g (the wrong combination) is expectedly much wider though it may contain 
some 'right' combinations due to accidental large fluctuation. 

Finally, we proceed to estimate the signal significance (S) using signal and back- 
ground events passing all the selection cuts mentioned above and use leading order cross- 
sections for signal as well as background processes. The effective cross-section for the 
signal events is 

aeff = a{pp titi*)5^ , 

where B is the total branching ratic^ for the i?-parity violating decays of the ti involving 
a 6-quark, namely ti — )■ bd, bs. It may be noted that the fraction of signal events passing 
the selection criteria increases with m^- although (T{pp — >■ titi ) falls rapidly. We choose a 
25 GeV wide window in the rn^^^^^^ distribution in which the signal is most abundant and 
has the right shape and determine the fraction of generated signal events {fsig) falling 
within this window. Finally the estimated number of signal events is obtained through 

Ns = fsig-0'eff--Cint , (3) 

^Note that limits from flavour changing neutral current processes constrain products of A" couplings [8]. 
The present study, though, is insensitive to how many channels share B. 
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Figure 6: The minimum value of a^jf (see text) required for excluding signal at 3cr for 
Lint = 100 fh~^ (solid line with open circles) and Lint = 300 fb~^ (solid line). The same 
required for observing a 5a excess over the background are shown for Lint = 100 fb~^ 
(dashed line with open circles) and Lint = 300 fb~^ (dashed line). For comparison, the 
solid line with open squares shows a{pp — )■ titi ). 



where Lint is the assumed value of integrated luminosity. Similarly, for each background 
process the number of events falling within the invariant mass window is determined and 
added up to obtain A'^^. The statistical significance may be estimated as 

S = . (4) 



Next, we use eqns.([3] & H]) to calculate the minimum value of effective signal cross- 
section {cr^ff) required to obtain S = 3 and S = 5 for a given luminosity (we use 

Lint = 100 fb~^, 300 fb~^). The values of titi production cross-section and cr^Iff have 
been listed in Table [1] and plotted in Fig. [6l It is evident that a large range of parameter 
space may be explored with Lint = 100 fb~^. 

Jet reconstruction and determination of invariant mass of a particle decaying into 
two or more jets become difficult at the LHC due to the presence of underlying events, hard 
radiation from the partons and large pt of the decaying particle. It has been claimed that a 
better invariant mass reconstruction is possible using more sophisticated jet reconstruction 
and analysis techniques and hence increase the signal significance [28]. As such techniques 
mature, they will help enlarge the scope of this analysis. 

In conclusion, we have explored the possibility of observing, at the LHC, a particle 
decaying into a pair of jets, of which one has a 6-quark as the originator. Concentrat- 
ing on the scalar top quark in a supersymmetric extension of the Standard Model, we 
demonstrate that simple and robust detector level observables may be used to explore a 
fairly wide range of the stop mass with integrated luminosities expected from the LHC. In 
addition, a fairly good estimate of the mass can also be obtained. And while our analysis 
has concentrated on a particular scenario, we have also demonstrated that it represents 
the most conservative choice amongst a class of such particles that appear naturally in a 
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Table 1: Table shows the minimum value of aeff (see text) required for different values 
of rrif^ (column 1) for exploring signal at 3(T for Lint = 100 fb~^ (column 3) and Lint = 
300 fb^^ (column 4)- The same required for observing a 5a excess over the background 
are shown for Lint = 100 fb~^ (column 5) and Lint = 300 fb~^ (column 6). Production 
cross- sections for the signal are shown in column 2. 



large variety of theories. Consequently, this analysis can be easily extended to such other 
scenarios allowing one to explore masses significantly larger than what this particular 
model allows us to. Aiding this assertion is the observation that, in many of such models 
such as those incorporating diquarks or excited quarks, the branching fraction into the 
relevant decay mode is overwhelmingly the dominant one. 
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